High Velocity Cloud Complex H: A Satellite of the Milky Way in a Retrograde 

Orbit? 



Felix J. Lockman 

National Radio Astronomy Observatory } P.O. Box 2, Green Bank, WV, 24944; 

jlockman@nrao. edu 

ABSTRACT 

Observations with the Green Bank Telescope of 21cm H I emission from the high- 
velocity cloud Complex H suggest that it is interacting with the Milky Way. A model 
in which the cloud is a satellite of the Galaxy in an inclined, retrograde circular orbit 
reproduces both the cloud's average velocity and its velocity gradient with latitude. The 
model places Complex H at R = 33 ± 9 kpc from the Galactic Center on a retrograde 
orbit inclined ~ 45° to the Galactic plane. At this location it has a mass in H I 
of > 6 x 10 6 Mq and dimensions of at least 10 x 5 kpc. Some of the diffuse H I 
associated with the cloud has apparently been decelerated by interaction with Galactic 
gas. Complex H has similarities to the dwarf irregular galaxy Leo A and to some 
compact high- velocity clouds, and has an internal structure nearly identical to parts of 
the Magellanic Stream, with a pressure P/k ~ 100 cm~ 3 K. 

Subject headings: Galaxy:halo — Galaxy: structure — ISMxlouds - Radio LinesdSM 



1. Introduction 

The current model for galaxy formation is one of hierarchical aggregation (Freeman & Bland- 
Bawthorn 2002), and in this context the problem of the nature of high- velocity B I clouds has 
attracted renewed attention (Blitz et al. 1999). The high- velocity B I cloud "Complex B", first 
discovered by Dieter (1971) and Bulsbosch (1971) and shown in Figure 1, covers more than 100 
square degrees at very low Galactic latitude. Its brighter components have a velocity with respect 
to the Local Standard of Rest of Vlsr ~ —200 km s _1 , which, in this direction, is at least 35 km 
s _1 more negative than is allowed for an object in circular Galactic rotation at any distance, and 
about 100 km s _1 more negative than the observed minimum velocity from R I in the Galactic disk. 
It lies in the Galactic plane, so the anomalous motion cannot be attributed to infall. 
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Wakker et al. (1998) show that Complex H must be more than 5 kpc from the Sun, which 
eliminates a possible association with the Local Bubble, and they dismiss several other "local" 
models on the grounds that the cloud does not seem to show strong interaction with Galactic gas. 
Morras, Bajaja, & Arnal (1998), however, note that Complex H overlaps in position with a large 
H I shell at permitted velocities, and propose that it created the shell as it penetrated the Galactic 
disk about 22 kpc from the Galactic Center. This interpretation has been challenged by Blitz et 
al. (1999), who cite a lack of ancillary evidence for so energetic an encounter and, instead, hold 
Complex H to be an important link in a chain of evidence which leads them to conclude that most 
high-velocity clouds are at intergalactic rather than interstellar distances. 

In fact, the basic conclusion that the velocity of Complex H cannot arise from circular rotation 
about the Galaxy is correct only for prograde orbits. An LSR velocity of —200 km s -1 at £ = 130° 
can be achieved by an object with Vg/R ~ —5 km s" 1 kpc -1 , where R is the distance from the 
Galactic center, Vg is the azimuthal velocity, and the negative sign indicates that the motion is 
against the sense of normal Galactic rotation. This Letter presents 21cm H I observations of 
Complex H made with the Green Bank Telescope which show that there is H I at velocities between 
the cloud and the disk, confirming previous suspicions that it is interacting with the Milky Way. 
The new data are consistent with a model in which this object is a satellite of the Milky Way 
in an inclined, retrograde orbit, whose outermost layers are currently being stripped away in its 
encounter with the Galaxy. 



2. Complex H and the Galactic Disk 

Figure 2 shows a cut through Complex H at t = 130? 5 made from the Leiden-Dwingeloo (LD) 
21cm H I survey data (Hartmann k, Burton 1997). The cloud divides into two distinct sections: a 
bright "core" which lies at b < 5°, and a fainter, more diffuse "tail" at b > 5°. Both the core and 
the tail move to more positive velocities away from the cloud apex near b ~ 4°, the tail reaching 
velocities allowed by Galactic rotation at b ~ 10°. This same general structure is seen at many 
longitudes and was apparent in the earliest observations (Dieter 1971; Hulsbosch 1975). There is 
also a hint of faint H I emission at velocities between that of the cloud core and the Galactic disk. 

Part of Complex H at 130?5 < £ < 132?5 and -2° < b < +8° was observed in the 21cm line 
with the Robert C. Byrd Green Bank Telescope (GBT) in August 2002 at 9' angular resolution 
with 3' spacing between spectra. Each spectrum covers —300 < Vlsr < +200 km s _1 at 1.25 km 
s _1 velocity resolution. The la noise in a channel is 0.16 K. The equipment and observing methods 
are identical to those described elsewhere (Lockman 2002). 

Figure 3 shows the GBT data at I = 130? 6. About 40% of the spectra through the core of the 
Complex have two line components, one broad and one narrow, with a median Av (FWHM) of 8 
and 28 km s _1 . In contrast, lines from the tail are almost uniformly broad with a median Av = 34 
km s _1 . The median width of all lines from the core is 22 km s _1 , while almost every spectrum 
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through the tail has Av > 22 km s . This suggests that the division between core and tail is not 
arbitrary, but is related to the nature of the Complex. Fig. 3 also shows diffuse H I emission filling 
much of the gap in velocity between the Complex and the Galactic disk, strong evidence for an 
interaction between the Complex and the Milky Way. The Complex is not nearly as isolated from 
the Galaxy as it seemed from earlier observations. 



3. A Model for Complex H 

In the model proposed here Complex H has a core which moves like a solid body on a circular 
orbit, and a tail of diffuse, decelerated gas resulting from its interaction with the Milky Way. The 
velocity of an object projected to the LSR can be written 



Vlsr = 



Vq Vn 

it XlQ 



cos(6) +y z sin(6), (1) 



in an R, 6, z polar coordinate system centered on the Galactic center, where 9 runs in the direction of 
Galactic rotation from the Sun-center line, and z is toward the North Galactic Pole. The standard 
Galactic constants, Ro = 8.5 kpc, Vq = 220 km s _1 are used throughout. If Complex H moves in a 
circular orbit (Vr = 0) about the Galactic center with a velocity of magnitude V c , at b = 0° it will 
have components Vq = V c cos(z) and V z = V c sin(i), where i is the inclination of the orbit to the 
Galactic plane: i = 0° is prograde in the Galactic plane, and i = 90° is toward the North Galactic 
Pole. 

The data constrain V z and Vq. At low latitudes cos(fe) ~ 1 and sin(fe) ~ b, so the gradient 
dVisR/db ~ V z . The observed slope in Vlsr of t ne core is about -3 km s _1 per degree of latitude 
giving V z = V c sin(i) = —170 km s -1 . Constraints on Vg/R come at b = 0°, where the second term 
of eq. 1 is zero. At £, b = 130?5 + 0° the core has V LS r = -195 km s" 1 implying that Vg/R = -4.3 
km s" 1 kpc -1 , or R = —0.23 V c cos(i) kpc. As the Galactic rotation curve is likely to be nearly 
flat to several hundred kpc from the Galactic center (Zaritsky 1999), setting V c = Vq = 220 km 
s _1 yields the orbital inclination i = 230° and the distance from the Galactic center R = 33 kpc. 
Complex H thus appears to be moving to greater longitude and to negative latitude at an angle 
w 45° with respect to the Galactic plane. The direction of motion is drawn in Fig. 1. 

These considerations are used to produce the simulation displayed in Figure 4. Here it is 
assumed that the Galaxy has a flat rotation curve with V c = 220 km s _1 and Complex H has two 
components, both with At; = 22 km s" 1 . The core follows the retrograde, inclined, orbit V c = 220 
km s _1 , i = 230°, at a distance R = 33 kpc from the Galactic center. The diffuse gas is slightly 
more than twice the size of the core and extends to higher latitude; its velocity varies linearly 
between the velocity of the core where they meet, and normal Galactic disk rotation, Vq = Vq and 
V z = 0, at its outer edge. This simulation is not intended to fit the structure of Complex H in 
detail, but simply to check the parameters of the proposed model. On this criterion the model does 
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well. The simulated core has the correct velocity and gradient, and the lagging diffuse gas produces 
a credible tail. 

Reasonable fits to the GBT data can be found for models in the range V z = —170 ±30 km s -1 , 
implying an approximate range of parameters i = 230 ± 10° and R = 33 ± 9 kpc. The orbit crosses 
the Galactic plane at 9 = 38 =L 5° from the Sun-center line in the direction of Galactic rotation at 
a distance from the Sun r = 27 ± 9 kpc. If the Galactic circular velocity declines beyond the Sun 
as is suggested by some evidence (Zaritsky 1999), a somewhat larger inclination would be derived, 
along with a smaller value of R. 

4. Implied Properties of Complex H 

At a distance from the Sun r = 27 kpc the total H I mass over Vlsr < —175 km s _1 (the shaded 
area of Fig. 1) is 6.4 x 10 6 M & . About half of this resides in the bright core at £ > 127°, b < +5° 
while about 10% is in the diffuse tail at b > 5°. There is certainly H I associated with the Complex 
at more positive velocities (e.g., that shown in the Fig. 1 contours), but the restrictions adopted 
here avoid possible confusion with unrelated material and provide a firm lower limit to the total 
mass. The size of the core in the £ and b directions is « 2 x 3 kpc, while the total extent of the 
cloud is at least 10 x 5 kpc. The dynamical mass needed to bind a spherical cloud of radius r c and 
line width Av is M^ yn « r c Av 2 /G, where G is the gravitational constant. The core of Complex H, 
with r c ~ 1 kpc and median Av = 22 km s _1 , has Md yn ~ 10 8 M Q , while the core H I mass is only 
3 x 10 6 M Q . Complex H is unlikely to be gravitationally bound by its H I alone. 

The two-phase structure of the H I lines from the core, also seen in other high-velocity clouds, 
suggests that it is in pressure equilibrium with a surrounding, low-density, medium (Cram & 
Giovanelli 1976; Ferrara & Field 1994). There is a group of H I features in the Complex at 
130° + 1° that have a characteristic angular size of 4' (Wakker & Schwarz 1991). At r = 27 kpc 
they would have a diameter of 30 pc and a typical (n) ss 2 cm~ 3 . Their temperature is likely 
50 K (Wakker, Vijfschaft, & Schwarz 1991) implying a pressure P/k = nT 100 cm~ 3 K. An 
identical pressure is derived from the GBT data, using the median core Av of 22 km s _1 to limit 
the temperature, and the median core Njjj of 6 x 10 19 cm~ 2 over a size of 2 kpc for an average 
volume density. 

Toward Complex H there is likely > 10 22 cm~ 2 of foreground H I. Studies of the region using 
the 2MASS survey have been complicated by strong differential reddening (0.6 < E(B-V) < 3), but 
no obvious stellar features attributable to Complex H are found in differenced, dereddened (J - K s , 
K s ) color-magnitude diagrams compared with nearby control fields (S. Majewski & M. Skrutskie, 
private communication). There is an IRAS point source toward Complex H tentatively identified 
as a possible young stellar object (Ivezic & Christodoulou 1997), and weak H a emission toward 
the brightest part of the Complex (Putman et al, in press), but the H a velocity is 30 km s _1 more 
positive than the velocity of the core and does not correspond to any feature in H I. 
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5. Summary Comments 

The model derived from the kinematics of Complex H is consistent with the overall structure 
of the object. The large-scale "head-tail" morphology seen in Fig. 1 (the trend continues to even 
more positive velocities than are shown) is the form one would expect from a cloud containing 
dense concentrations moving on an inclined, retrograde orbit to higher t and lower b, shedding 
gas through interaction with the Milky Way. A number of high-velocity clouds believed to be 
interacting with the Milky Way also have this structure, where, like Complex H, the velocity of the 
diffuse gas "downstream" in the tail is closer to Galactic rotation than the gas in the head (Pietz 
et al 1996; Brims et al. 2000). The physical issues involved in the interaction of a high- velocity 
cloud with an external medium have stimulated considerable theoretical interest (Benjamin 1999); 
magnetic fields may be important in preserving the integrity of the main cloud (Santillan et al 
1999; Konz, Bruins, & Birk 2002; Quilis & Moore 2001). 

In its small-scale structure Complex H resembles parts of the Magellanic Stream. Though its 
origin is tidal, the Stream shows evidence for interaction with an external medium (Putman et al 
1998; Konz, Brims, & Birk 2002). It contains H I in two phases, and clouds with a size of 15-30 
pc, (n) ~ 2 — 6 cm~ 3 , and Av ~ 5 km s _1 , virtually identical to the knots in Complex H (Wakker, 
Oosterloo, & Putman 2002). In its large-scale characteristics, Complex H resembles the dwarf 
irregular galaxy Leo A. At a distance of 0.8 Mpc it has about the same H I mass as the Complex 
and the same linear size as the core; the H I lines from the two objects have a similar width and 
two-component structure (Young & Lo 1996; Schulte-Ladbeck et al 2002). Most dwarf irregular 
galaxies, however, have considerably more H I mass than Complex H (Mateo 1998). The Complex 
resembles as well some "compact" high- velocity H I clouds which, if they are at distances of several 
hundred kpc, have a similar size and H I mass and also have two H I phases (Burton, Braun, & 
Chengalur 2001; de Heij, Braun, & Burton 2002). 

Complex H may be interacting with the Galactic halo or some part of the disk. Its orbital 
period for V c = 220 km s~ 1 is ~ 10 9 yr over which time it traverses heights of ±25 kpc from the 
Galactic plane. There is convincing evidence for a gaseous halo that extends to at least R ~ 50 
kpc toward the Magellanic Stream, though little is known of it in detail (Savage 1995; Kalberla 
& Kerp 1998; Konz, Brims, & Birk 2002; Savage et al 2003). The normal H I disk emission at 
t s=s 130° goes to Vlsr ~ — 110 km s" 1 which, for a flat rotation curve, is R = 25 kpc, similar to the 
derived distance of the Complex. The Galaxy may also have an extended disk of ionized gas like 
some other systems (Bland-Hawthorn, Freeman, & Quinn 1997). The existing GBT H I data do 
not cover enough area to be useful in evaluating the suggestion of Morras, Bajaja, & Arnal (1998) 
that Complex H has actually pierced the H I disk, but further observations are planned. 

Galaxies exist in a rich environment of companions, interaction and distortion (Sancisi 1999). 
Retrograde motions are observed in H I components of many systems (Bertola & Corsini 1998), and 
in Milky Way globular clusters (Freeman & Bland-Hawthorn 2002). A retrograde orbit at large R 
tends to be less destructive to a satellite galaxy than a prograde orbit (Toomre & Toomre 1972; 
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Hernquist & Quinn 1988), which may account for the persistence of the object. But Complex 
H may not persist much longer. If a complete accounting of the mass of diffuse, presumably 
decelerated, gas associated with the Complex could be made, it might well show that the coherent 
core is only a fraction of the total H I, implying that Complex H is now in the process of being 
destroyed. 

Finally, the guiding assumption that Complex H is in a purely circular orbit (or is just now at 
perigalacticon) is clearly too naive — there is likely some ellipticity and a radial velocity component. 
Additional data will supply extra constraints on Vq. But the simple model presented here, motivated 
by the new evidence for interaction with the Galaxy, suggests that the velocity of Complex H and 
its gradient with latitude carry significant information about its space motion, and are not entirely 
contingent. 

I thank R.A. Benjamin, W.B. Burton, S. Majewski, E.M. Murphy, T.M. Tripp, and the referee, 
H. van Woerden, for useful discussions and comments on the manuscript. 
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Fig. 1. — Complex H in the Leiden-Dwingeloo 21cm H I survey data (Hartmann & Burton 1997) 
at 35' angular resolution. The shading shows H I at Vlsr < —175 km s _1 , which includes the 
brightest H I in the Complex. The contours trace the H I emission at more positive velocities 
between —175 and —140 km s -1 , closer to that allowed by normal Galactic rotation. The vertical 
lines near I = 130° show the locations of the velocity-latitude cuts of Figs. 2 and 3, which were 
used to derive the model for the cloud's motion in the direction shown by the arrow. In this model, 
most of the contoured material would have been decelerated by interaction with the Milky Way. 
The grey scale is in units of 10 19 cm -2 ; coutours are equally spaced every 2 x 10 19 cm -2 . 
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Fig. 2. — A velocity-latitude cut (a "slit spectrum") through Complex H at I = 130? 5 from the LD 
survey data showing the "core" at b < 5° and the "tail" at b > 5°. The slanting line through the 
core follows Vlsr = —195 — 170 sin(6) km s _1 . The slope of the line gives the vertical component 
of orbital motion, V z , while its value at b = 0° constrains Vg/R. Contours are in factors of two 
from the peak; the gray scale shows emission at T& < 3 K. 
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Fig. 3. — A velocity-latitude cut through Complex H at I = 130?6 made from GBT data at 9' 
angular resolution. The two panels show identical data; the upper panel is clipped at 5 K, the 
lower at 1 K. The "tail" at b > 5° is almost entirely broad- line gas while the bright "core" at b < 5° 
contains many knots of narrower lines. There is diffuse H I emission at velocities between Complex 
H and the Galactic disk suggesting that some material from the Complex has been decelerated by 
interaction with the Milky Way. 



- 12 - 





o 
m 



o 
o 




o 

LO 



o 
o 



" o 
in 

i , , , , i , , , , i , , , , i , , , , d ^ 

01 S 5- 01- 

8pn}l}D"| 0ipD|D3 



Fig. 4. — Simulation of Complex H as a cloud in an inclined retrograde orbit surrounded by a 
diffuse component whose velocity changes linearly between the retrograde core and normal prograde 
Galactic rotation. The simulation is intended to match the observations of Figs. 2 and 3. The 
velocity tilt of the core results from the projection of its vertical motion. The diagonal line is the 
same one drawn through the data in Figure 2. 



